Ferromagnetic spin correlations in a few-fermion system 
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We study the spin correlations of a few fermions in a quasi one-dimensional trap. Exact di- 
agonalization and complementary Quantum Monte Carlo calculations demonstrate that repulsive 
interactions between the two species drives ferromagnetic correlations. We show how the ejection 
probability of an atom provides an experimental probe of the magnetic spin correlations. With more 
than five atoms trapped, the system approaches the itinerant Stoner Hmit. Finally, studying the 
energy levels of the excited states reveals that losses to Feshbach molecules are suppressed by the 
discretization of energy levels when fewer than seven atoms are trapped. 

PACS numbers: 67.85.Lm, 03.65. Ge, 03.65.Xp 



Recent experimental advances allow investigators to 
confine up to twenty atoms in a single trap and ad- 
dress their quantum state [3 [5] . This precise control en- 
abled the Heidelberg group to confine two distinguishable 
fermions in a one-dimensional well and tune the interac- 
tion strength to induce fermionization [3], allowing the 
exact analytical confirmation of the fundamental physics 
of short range repulsion |1HT5]. Here we demonstrate 
that just such repulsive interactions acting between a few 
fermions allows us to construct a Hamiltonian analogous 
to the Stoner model |14j and offers experimentalists an 
opportunity to observe the emergent ferromagnetic cor- 
relations without losses to Feshbach molecules. 

The itinerant ferromagnet predicted by the Stoner 
Hamiltonian has never been cleanly realized and stud- 
ied in the solid state. However, it was proposed [TSl - 
[TT] that a ferromagnetic phase could be created with 
a fermionic cold atom gas. An experiment on a ^Li 
gas by the MIT group displayed signatures compatible 
with ferromagnetism [TSlEn] . but the observations were 
later explained by a loss mechanism pTl - 125] . To circum- 
vent losses it has been suggested that magnetic corre- 
lations could be explored in systems with a mass im- 
balance [23], two-dimensional geometry [5S], spin spi- 
rals [35], or flux lattices [17] • Here we demonstrate how 
a quasi one-dimensional system containing only a few 
fermions could avoid the problems associated with losses 
and display ferromagnetic correlations. Our main result 
is shown in Fig. IT] the discretization of energy levels in 
the few fermion system means that losses to Feshbach 
molecules are restricted to a narrow range of interaction 
strengths, allowing a tunneling measurement of the ejec- 
tion of an atom to expose the underlying ferromagnetic 
correlations. 

In this paper we first describe the proposed experimen- 
tal setup, then demonstrate the emergence of ferromag- 
netic correlations that we study through a tunneling pro- 
cess. Finally, we demonstrate that the formation of the 
competing dimer state is inhibited by the discretization 
of the energy levels in the harmonic confining potential. 



FIG. 1. (Color online) Left: The green solid line is the trap- 
ping potential V(z) that is lowered to the green dashed line 
to allow an atom to escape. The energy levels for up-spin 
(red) and down-spin (blue) atoms within the potential are 
shown in brown, with the putative escape of the highest en- 
ergy level up-spin atom. Right: The probability of ejection of 
a down-spin atom into a state fully polarized along the quan- 
tization axis. The gray region is excluded due to atom losses 
for differing initial numbers of atoms. 



EXPERIMENTAL SETUP 

A fermionic gas of two hyperfine states with pseudospin 
a € {t, i} is tightly confined as shown in Fig. flj We seek 
to solve the full Hamiltonian in absolute coordinates H ~ 
~n^V'^l2m + muj]_{x^ + y'^)/2 + mLj'^^z'^/2 + F(ri - ra), 
with m the atomic mass. The confining potential is ax- 
ially symmetric with trap frequencies uj± = Tywy and 
T] ^ 10 [21 [3], and we define the harmonic oscillator 
lengths Qi = y/2K/rnuji. Only a single transverse mode 
is occupied, constraining the atoms into the quasi-lD 
regime. In this limit we re-parameterize the interspecies 
potential V{r) = —U&{R— |r|) first into the s-wave scat- 
tering length a3D — R[l — ta,n{x) / x] with x = RvmU/h, 
and then as a one-dimensional pseudopotential [28] g = 
h'^azT>/ma^{a^--Ca^D) with C = C,{l/2) « 1.46. A con- 
finement induced resonance emerges at osd — a\^IC We 
verified that the results tend to the contact limit below 
R = 0.2a||. In the limit aj|| — > 0, we recover the Stoner 
Hamiltonian H — —t?\l'\l2ra + tocjjJz^/2 -|- g8[z\ — Z2). 

To probe the quantum state we apply a magnetic field 
gradient to tilt the external potential (see Fig. fl]) and al- 
low one atom to escape. The remaining atoms can then 



be characterized by the single atom addressabihty [2]. 
We denote the number of trapped spins N^ and N^ so 
that investigators can directly measure the spin in the 
quantization direction S^ = {Ni^ — -^^)/2 and the to- 
tal number of atoms N^ + N^ in the final state. How- 
ever, our measure of ferromagnetic correlations, the spin 
S = (I2ni4it cl;)''^'(cnT Cni)"^) Can bc aligned in any 
direction, where a denotes the vector of Pauli-spin ma- 
trices and Cjia- is the annihilation operator of an atom of 
spin a from harmonic oscillator state n. The Hamiltonian 
is SU(2) invariant so the spin quantum number defined 
through s{s + 1) = (S^) is a good quantum number, al- 
lowing us to define the quantum state \s,N^,Ni). With 
two atoms a polarized s = 1 state can be generated not 
only from the S^ — I state, denoted 1 1, 2, 0) , but also from 
the Sz — state denoted |1,1,1), corresponding to the 
prototypal triplet states |tt) and (|ti) + |it>)/%/2- The 
unpolarized Sz = state, denoted |0, 1,1), corresponds 
to the singlet state (|ti) - |;t))/\/2. 



ENERGY OF STATES 

Our main tool to study the system is exact diagonal- 
ization. We build a one-atom basis from the Gaussian 
orbitals 4>n^^n .n^ix,y,z) of the harmonic trapping po- 
tential. We retain all orbitals that satisfy (nx,%) < 2 
and Uz < 20. We construct the Slater determinants for 
this basis and select the 10,000 determinants with low- 
est non-interacting energy, guaranteeing convergence of 
the energy of the open channel within O.OObHujH at the 
confinement- induced resonance. This is much smaller 
than the energy scale of magnetization, N^NiHui^^ and 
less than the energy difference from the true itiner- 
ant state shown in Fig. l2jf). We next calculate the 
interaction matrix elements numerically, construct the 
Hamiltonian matrix for a particular g using the Slater- 
Condon rules, and diagonalize it to obtain the eigenstates 

{^m(ff)}- 

Exact diagonalization is restricted to systems of fewer 
than five atoms. However, the experimental setup can 
contain up to twenty atoms so to analyze the general 
many-body case we use the QMC code CASINO ^51150] . 
The approach is a refinement of that used in previous 
studies of itinerant ferromagnetism |THl [5TH55] . We use 
a trial wave function ip — FD that is a product of 
a Slater determinant, D that takes full account of the 
fermion statistics and a Jastrow factor F that includes 
further interparticle correlations. We use the lowest- 
order constrained variational method [551 [57] that is 
common in nuclear physics and has also been used to 
study cold atom gases including ferromagnets [SH |35l 
[Mll^- This method solves the Hamihonian [-d^/dr^ + 
mV{r)]rf{r) — k'^rf{r). For low energy s-wave scat- 
tering this gives /(r) w 1 — a^u/'r that has a node at 
the scattering length osd, and saturates at large dis- 
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FIG. 2. (Color online) (a-d) The energy of the (two-five)- 
atom states calculated with exact diagonalization. The red, 
magenta, and blue lines highlight open channels with § signi- 
fying identical states with larger S^, and gray lines indicating 
all other states. The green line shows the final bound molecu- 
lar state that crosses the low-spin open channel, and the green 
dashed line the bound state with the COM motion excited 
into the nx,y = 2 orbital examined in Ref. [13]. The points 
show QMC results whose uncertainty is the point size, (e) 
Density profile of the state |l/2, 2, 1) for two majority species 
atoms with a pinned minority atom, (f) Difference in energy 
from the polaron in an infinite system ,43j with number of 
trapped atoms, (g) Pair correlation function at the confine- 
ment induced resonance around a down-spin pinned at z = 
normalized by the non-interacting density profile. 



tances. This solution is embedded into a Jastrow factor 
F = Y\i j f{\^i—i^j\) where r^ is the position of the ith up- 
spin and Tj the jth down-spin to guarantee occupation of 
the upper branch |34l 135] . In the quasi one-dimensional 
setting the transverse confinement does not allow occu- 
pation of the bound state since its maximum binding 
energy ^ — 1.5fi.cij|| is much less than the energy scale 
of occupation of transverse modes ~ 20huj\i . This was 
further confirmed by studying the exact diagonalization 
states, where deep in the Super- Tonks regime occupation 
of the higher transverse modes to be ^ 10~^, resulting in 
the strong agreement between exact diagonalization re- 
sults and QMC demonstrated in Fig.l2ja-d). Meanwhile, 
to calculate the binding energy of the molecule a.t g > 



and the ground open channel state at g < we use a 
Jastrow factor F = e"^ , where J includes the polynomial 
expansion in atom-atom separation proposed in Ref. [40] 
with eight variational parameters that we optimize. 

With the exact diagonalization and QMC formalism 
in place, in Fig. [2ja-d) we compare the ground state en- 
ergy predicted by both exact diagonalization and QMC, 
and also the two atom exact analytical solution [H [6l [7] . 
There is strong agreement at all interaction strengths. 
The underlying attractive potential means that exact di- 
agonalization also delivers the multitude of molecular 
bound states and repeated bands incremented by Huj^^ 
corresponding to center-of-mass (COM) motion. We 
now focus on the two atom system Fig. ^a) where 
the open channel |0,1,1) can be adiabatically tracked 
across the bound states. In the non-interacting regime, 
— ScLiyaii/g — >■ —CXI, the lower spin state |0, 1,1) has the 
lowest energy. At the confinement induced resonance the 
spin states cross [21 E], and in the Super- Tonks regime, 
—Hu!\\a\\/g —7' oo, the s = 1 states have lower energy. 

In the three-body system in Fig. l2Fb) three open chan- 
nel states are possible: the low spin 1 1/2, 2,1) and the 
high spin states |3/2,3,0) and |3/2,2, 1). Similarly to 
the two-body system, at weak interactions the s = 1/2 
state has the lower energy, the bands cross at the confine- 
ment induced resonance, and in the Super- Tonks regime 
the s = 3/2 states are favorable in good agreement with 
existing literature [HI [H]. In Fig. ^c) we also studied 
the four atom case where three values for the spin are 
available: s E {0,1,2}. The three bands cross at the 
confinement induced resonance meaning that any poten- 
tial onset of ferromagnetism would be abrupt, as occurs 
in the infinite body case [HJ HI] . With five atoms the 
molecular bands become more prevalent, and we will later 
demonstrate how they allow losses into bound molecules. 

In the Super- Tonks regime the high spin state is ener- 
getically favorable so the gas would enter the magnetic 
phase if it were not blocked by spin conservation. How- 
ever, the spatial distribution of the atoms betrays the 
underlying magnetic correlations. We focus here on the 
three atom state |1/2,2,1), and pin the minority down- 
spin atom at z — — 1.5aj|. In the non-interacting case 
the up-spin atom density is concentrated around the trap 
center irrespective of the down-spin position. At the con- 
finement induced resonance the up-spin density is driven 
to zero at the down spin pinning position, whereas in the 
Super- Tonks regime, the up-spin atoms are forced away 
from the down-spin forming a separate magnetic domain. 

Now that we have observed the emergence of mag- 
netic correlations we study the energy of a single down- 
spin in a trap with iV-|- majority spin atoms to assess 
the consequences of system size and whether the system 
can be regarded as a model for the Stoner Hamiltonian. 
We compare our system to the analytically solvable po- 
laron limit [43) of a single down spin in a sea of up-spin 
atoms. We first study the energy of a polaron in Fig.[2j 



With N^ — 1 the exact diagonalization agrees with the 
Busch analytical result [7] , showing less repulsion energy 
than the infinite body case. On inserting more major- 
ity spin atoms the energy quickly tends to the infinite 
sized limit [13], being within < 1% at all interaction 
strengths with A'^tot > 5. Secondly, in Fig. [2]jg) we study 
the pair correlation function. With more majority spin 
atoms the pair correlation function quickly tends to the 
infinite body limit [33], with the correct correlation hole 
and first Friedel oscillation observed for A'tot > 5. Both 
pieces of evidence indicate that systems with Atot > 5 
can be regarded as reasonable representations of the itin- 
erant Stoner Hamiltonian. 



TUNNELING STATISTICS 

Although the magnetic phase is energetically favorable 
in the Super- Tonks regime its formation is prohibited by 
spin conservation. In Fig. IT] we therefore consider tilt- 
ing the trap to allow one atom to escape. This allows 
the system to tunnel into the magnetic ground state con- 
taining one fewer atom. We calculate the tunneling rate 
using Fermi's Golden rule. The tunneling rate F exhibits 
an exponential dependence on the escape energy, such 
that we need only consider tunneling from the highest 
occupied orbital with maximal energy -Eesc- We now 
consider a general intermediate interaction strength and 
calculate the probability of forming a particular state i. 
Pi = rj/^.Fj. This tunnehng probability calculated 
from the exact diagonalization data is shown in Fig. IT] 
We focus on two contrasting starting states: the polaron 
limit with multiple up-spin atoms and a single down-spin 
atom, and the unpolarized limit with equal numbers of 
the two species. 

Polaron: At zero interactions the highest energy ma- 
jority spin atom is expelled leading to zero probability of 
ejecting the minority spin atom, whereas at the confine- 
ment induced resonance all atoms have an equal proba- 
bility of expulsion. In the Super- Tonks regime starting 
with Atot atoms, the system will tunnel most rapidly into 
the state with lowest energy - the fully polarized state 
with s — (Atot — l)/2- There are two quantum states 
available: with Sz = (A^tot — l)/2 formed by the ejection 
of the down-spin atom and Sz = (A^tot ^ 2)/2 formed 
by the ejection of an up-spin atom. Starting with three 
atoms tunneling into a two atom state, these are simply 
the triplet states |tt) and (|tl) + ||t))/\/2. These two 
possibilities occur with equal probability, giving a plateau 
probability of 1/2 for ejecting a minority spin atom in the 
Super- Tonks regime irrespective of the initial number of 
atoms. The probability curves in Fig. [l] become increas- 
ingly sharp with a higher number of atoms because of the 
larger energy exchange over the same range of interaction 
strengths. 

Unpolarized limit: An alternative experimental proto- 
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FIG. 3. (Color online) (a) The energy bands, following Fig. [2] 
conventions, (b) The map of crossing events for the riband- 
th excited molecular state with NfNi. Blue denotes a single 
crossing, and red a double crossing, (c) The minimum inter- 
action strength to avoid the band crossing region. 



col is to start from an S^ — state and eject half of the 
atoms. In the Super- Tonks regime the system will pass 
through a series of partially polarized states with a final 
probability of being fully polarized along S^ of 2i^^t°'/2. 
As seen in Fig. IT] starting with four atoms this is 1/2, 
whereas with six starting atoms it is 1/4, providing an 
ideal test of the final state statistics. 

The tunneling method is sensitive to S^ but not s so 
does not provide a full diagnosis of the final quantum 
state. This is exemplified when starting from the po- 
laron state in the Super- Tonks regime where the ejection 
probability of a minority spin is 1/2 rather than unity. 
To distinguish the 1 1,1,1) state from the other possi- 
ble ^z = state, |0, 1, 1), one could ramp the interaction 
strength into the non-interacting regime and measure the 
energy of the state through a second tunneling measure- 
ment [HIS]. Should the |1, 1, 1) state be dominant the en- 
ergy will be independent of interaction strength, whereas 
if |0, 1, 1) dominates, the energy will fall. 



LOSS MECHANISM 

The search for itinerant ferromagnetism in a cold atom 
gas has been plagued by a competing loss process piH^ . 
Several models for loss have been put forward including 
two and three-body models [5T1 1^ HH US], and losses 
to states excited with transverse COM motion [131 ES] . 
Both mechanisms require the crossing of bound molecu- 
lar bands to obey energy conservation. Our Hamiltonian 
only displays avoided crossings between states with the 
same COM quantum number. However, states with dif- 
ferent COM motion could have avoided crossings due to 
unforeseen perturbations such as an anharmonic poten- 
tial [131 123]. To guarantee a loss-free experiment, we 
fence off the region in which the open channel is crossed 
or anti-crossed by any other state. This pessimistic ap- 
proach will be robust to unforeseen perturbations that 
may alter the crossings but will not significantly alter 
the positions of the bands. We first focus on the three- 
body system where we use Fig.[2]jb) to define a loss region 
as where the desired open channel |l/2, 2, 1) crosses the 



molecular bound states. 

The ground molecular bound state labeled (i) is lower 
than the entire open channel |l/2, 2, 1) so its formation is 
prohibited by energy conservation. The molecular bound 
state can be excited with COM motion, giving rise to in- 
creasingly populous families of curves. The curves (ii) are 
the final set that do not cross the state |l/2, 2, 1), whereas 
the family of curves (iii) do cross the state 1 1/2, 2,1). 
Crossings allow losses to molecules that conserve energy, 
so borders a region of molecule formation. Further cross- 
ings from more excited molecular bands occur up to 
the confinement induced resonance, prohibiting experi- 
mentalists from looking for magnetic correlations within 
< —huJi\ai\/g < 0.24. The region contains the molecular 
bound state with COM motion excited into the second 
transverse mode that is a significant cause of loss in an 
anharmonic potential. [131 Si] 

We note that it is possible to adiabatically transit 
across the region of band crossing. Investigators can per- 
form experiments on the |l/2, 2,1) state either side of 
the shaded region in Fig. IT] but not within it. A similar 
analysis of a system with four atoms reveals that losses 
would block the region < —HcoHaH/g < 0.36, and with 
five atoms the range < —hw^aw/g < 0.42. 

Exact diagonalization cannot accurately address larger 
systems. We therefore turn to QMC for the open channel 
and the variational QMC for the molecular band. In the 
Super- Tonks regime the energy difference between unpo- 
larized and polarized states is hujuN^N^, so we categorize 
states by N^N^. We focus on the state that bounds the 
loss region, with the molecule having no COM motion, 
and other atoms in higher energy orbitals compatible 
with the correct non-interacting energy. When we had 
N^N^ = 2 in Fig. [2jb) the excited molecule bands cross 
the upper branch only once, whereas with N^N^ = 8 in 
Fig. IsFa) a molecular band, highlighted in green, crosses 
the open channel twice. In Fig.lslb) we show whether the 
nband-th family of excited molecular bands crosses the 
open channel once or twice. The double crossings first 
emerge at N^N^ = 7 and become ubiquitous as N^N^ 
rises. This leads to a proliferation in the total number of 
band crossings, and as shown in Fig. Islc) a dramatic rise 
in the minimum interaction strength — fia||a;||/g required 
to avoid band crossings. 



DISCUSSION 

A quasi one-dimensional system containing a few 
fermionic atoms poses an opportunity to explore ferro- 
magnetic correlations. Discretization of the energy lev- 
els means that the system offers the stabilization of a 
ferromagnetic phase without losses. We have calculated 
the energy structure and studied the experimentally ob- 
servable ejection probabilities. We have demonstrated 
that both the polaron energy and pair correlation func- 



tion tend to the itinerant limit when TVtot > 5, whereas 
molecule losses restrict the experimental observation of 
magnetic correlations to -/Vtot < 6, meaning that the sys- 
tem with iVtot S {5, 6} could present an opportunity to 
observe magnetic correlations driven by the Stoner mech- 
anism without losses. 
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